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Abstract 

We propose simple heuristics for the assembly line worker assignment and 
balancing problem. This problem typically occurs in assembly lines in sheltered 
work centers for the disabled. Different from the classical simple assembly line 
balancing problem, the task execution times vary according to the assigned worker. 
We develop a constructive heuristic framework based on task and worker priority 
rules defining the order in which the tasks and workers should be assigned to the 
workstations. Wc present a number of such rules and compare their performance 
across three possible uses: as a stand-alone method, as an initial solution generator 
for meta-heuristics, and as a decoder for a hybrid genetic algorithm. Our results 
show that the heuristics are fast, they obtain good results as a stand-alone method 
and are efficient when used as a initial solution generator or as a solution decoder 
within more elaborate approaches. 

1 Introduction 

In assembly lines, products are assembled by means of the successive execution of 
tasks in workstations. Each task has an execution time and precedence relationships 
with other tasks. The simple assembly line balancing problem (SALBP) concerns the 
decision of allocating tasks to the workstations while respecting the partial ordering of 
the tasks. Let c be the cycle time of the line, i.e., the time spent at one of the worksta- 
tions with the heaviest workload and m be the number of workstations. The problem 
is known as SALBP- 1 when the goal is to minimize m given a maximum allowed c, 
SALBP-2 when the goal is to minimize c given a fixed m, SALBP-E when one wants 
to minimize m ■ c and SALBP-F when the goal is to decide whether a feasible solution 
does exist for given m and c. For more details on SALBP and many of its variants, we 
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refer the reader to |Baybars| ^T9S6[, |Scholl| ( [1999) [Scholl and Becker] ( [2006) 1, [Boysen 
etal[ ( [2007) 1, and jBoysen et al| ( |2008{ 



The assembly line worker assignment and balancing problem (ALWABP) is an ex- 
tension of the SALBP in which task execution times are worker-dependent. This prob- 
lem typically occurs when balancing assembly lines with disabled workers, because 
a given worker might be efficient on a certain subset of the tasks while being ineffi- 
cient on (or even unable of executing) other tasks. Analagous to the SALBP, variants 
ALWABP- 1, -2, -E and -F can also be defined. 

The ALWABP has been introduced in the literature by jMiralles et al ( 2007| l by 
means of a case study in an assembly line of a Spanish sheltered work center for the 
disabled where a fixed number of workers should be present in an assembly line whose 
production rate should be maximized (ALWABP-2). The same authors have later de- 
veloped a branch-and-bound algorithm for the problem, enabling the solution of small- 
sized instances ( [Miralles et al||2008 ). Because of the problem complexity and the need 
to solve larger instances, the Uterature has since then shifted its efforts to heuristic 
methods. [Chaves et al ( 2009] ) have proposed an elaborate clustering search algorithm 
that obtained good solutions in reasonable computation times for instances of up to 19 
workers. Their algorithm relies on a clustering procedure to group, analyze and im- 
prove the solutions generated by a meta-heuristic approach. In the search of a simpler 
method, Moreira and Costa ( 20091 have proposed a minimalist tabu search algorithm 
based on the use of objective function penalties associated to infeasible solutions to 
improve the algorithms ability to effectively explore the search space. More recently, 
Blum and Miralles (201 U have developed a beam search algorithm which has been 



successful in improving the best known solutions obtained in the previously mentioned 
studies. 

In this study, we develop and test constructive heuristics for the ALWABP-2. The 
heuristics generalize station-oriented assignment procedures for the SALBP and ap- 



ply them in search procedures for finding the minimum cycle time ( Scholl and VoB 
1996| l. Station-oriented assignment procedures sequentially process the stations and 



assign tasks to the current station in an order defined by a priority rule until it is max- 
imally loaded. The extension of such procedures to the ALWABP-2 is not immediate, 
since the priority rules are often based on the task execution times, a parameter that in 
the ALWABP depends on the worker assigned to each workstation which is itself an 
optimization decision. 

The idea of our method is to use task and worker priority rules to define which 
worker and which set of tasks will be assigned to each workstation. We define 16 task 
priority rules and three worker priority rules. The heuristic efficiency is compared by 
means of computational tests to evaluate its performance both as a stand-alone method 
as well as its ability to improve the convergence of more sophisticated methods such as 
the aforementioned clustering search and tabu search algorithms ( [Chaves et al| [2009 [ 
[Moreira and Costa[|2009) ). Moreover, the proposed heuristic is also used as a solution 
decoder within a hybrid genetic algorithm that optimizes explicit priorities for each 
task- worker pair. 

The remainder of this paper is organized as follows. Next, we present a formal 
definition of the ALWABP and a mixed-integer linear model for the problem. Then, the 
proposed heuristic is detailed in Section[3[and its use within a hybrid genetic algorithm 
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is described in Section]?] A computational study is presented in Section |5] and this 
paper ends with some conclusions in Section|6] 



2 Formal problem definition and mathematical model 

In this section we present a formal definition of the ALWABP-2. As mentioned before, 
we denote by c the cycle time and by m the number of workstations. Let be the set of 
tasks to be allocated. The precedence constraints are given by a directed acyclic graph 
G — {N,E) over this set of tasks, where each edge ij € E indicates that task / is an 
immediate predecessor of task j. We also define an extended version of the precedence 
graph known as the transitive closure of G, G* = {N,E*), in which ij £ E* if there is a 
path from / to /' in G. In addition we will use the following notation in the remainder 
of this paper: 

5 set of workstations; 

W set of workers, |W| — \S\; 

ti e N time of task / (SALBP); 

tn-i G N U {00} time of task / when executed by worker w (ALWABP); 

Iw = {i & N \ t„i = 00} set of tasks infeasible for worker w; 

Pi = {i\ ji £ E} set of immediate predecessors of task /; 

P* = {j \ ji EE*} set of all predecessors of task /; 

Fi — {j \ ij G E} set of immediate successors of task /; 

F* ~ {j \ ij EE*} set of all successors of task /. 
The goal of the problem is, given a fixed number of workstations, m, to find an 
assignment of tasks and workers to the workstations minimizing c and such that each 
task i E N is assigned to a single station s E S and the precedence relationships are 
respected, i.e., a task / E Pj can only be assigned to the same workstation to which 
task j has been assigned or to workstations preceding it. We consider an equal number 
of workstations and workers and, therefore, each worker w eW must be assigned to 
a single station s E S. Likewise, each workstation s E S must receive a single worker 
w EW. Workers cannot be assigned independently from the task assignments, since 
the total load of a workstation depends both on the tasks and the worker assigned to it. 

A mathematical model has been proposed by |Miralles et a l (2008). The authors 
define binary variables x^m, equal to 1 only if task ; is assigned to worker w at worksta- 
tion s, and binary variables ysw, equal to 1 only if worker w is assigned to workstation 
s. The model can be written as below: 

Minimize c (1) 

subject to = ^ V/€A^, (2) 

weWseS 

£}'™. = i ywEW, (3) 

seS 

£ Xvv. = 1 e 5, (4) 

wew 

£ ^xsni < L L -'^^"v V/, jEN\iEPj, (5) 
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'v"-*^™' Vi e 5, (6) 

wew ieN 

Y -^.»«- < l^b^-w ewyse S, (7) 

ieN 

x,„i = Vvv e w, Vi e 5, V/ e /„. , (8) 

y,„,e{0,l} Vie5,VweW, (9) 

xswi e {0, 1 } e 5, Vw e w, V/ e A?. (10) 

Model ([T]i-([TO| focuses on minimizing the cycle time for a given number of work- 
stations. Constraints (j2]) guarantee that each task is executed, and that it is done by 
a single worker, at a single workstation. Constraints ([3]) and (|4]) establish a bijection 
between workers and workstations at a feasible solution, i.e., every worker is assigned 
to a single workstation and vice versa. Constraints (jSj) define the precedence relations. 
Constraints (j6]l establish that the cycle time is the sum of the execution times of the 
tasks at the most charged workstation, by forcing c to be larger than or equal to the 
load assigned to each station. Finally, constraints (|7]l indicate that a task can only be 
assigned to a worker in a given workstation if that worker is also assigned to the work- 
station. 



The difficulty in solving this model has been reported in the literature (Miralles 
|et al||2008 ^) and motivates the search for efficient heuristic methods like the constructive 
heuristics presented in the following section. 



3 Constructive heuristic methods 

In this section we propose a constructive heuristic framework for the ALWABP which 



is based on heuristics for the SALBP presented by SchoU and VoB ( 1996]l. In order to 



simplify the presentation, we describe the original heuristic in Section [3T| and present 



the new developments in Section 3.2 



3.1 Constructive heuristics for the SALBP 

The strategy used by |Scholl and VoB| ( |1996| l to solve the S ALBP-2 relies on solving 
instances of SALBP-1 for different cycle times. A given cycle time c can be considered 
an upper bound for the S ALBP-2 if the solution of the SALBP-1 needs, at most, the 
desired number of workstations m. 

The solution of the SALBP-1 instances can rely upon exact or approximate meth- 
ods. Since the SALBP-1 is also an NP-hard problem, heuristic methods are often used. 
These methods are usually based on priority rules that order the tasks according to a 
given criterion and assign them to the workstations accordingly. An important notion is 
that of an available task, which is an unassigned task whose predecessors have already 
been assigned. Let c be the desired cycle time and f (5^) the station time of workstation 
k. A station-oriented procedure to solve the SALBP-1 can be easily described in three 
steps: 



Scholl and VoB ( 1996 1 mention the following selection of successful priority rules. 



defining different orders with which the tasks should be considered for assignment: 
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Algorithm 1 Station-oriented assignment procedure for the SALBP-1 
Require: maximum cycle time c 
1: Start at an empty workstation A; = 0, /(5j.) = 0. Set the unassigned tasks to U = N. 
2: Assign to the current workstation the available task / with highest priority such that + 

t{Sk) <cand setU = U\{i}. 
3: If no more tasks are available, return k (end). Otherwise, if there are available tasks, but 
none can be assigned to the current workstation k, set k = k+ 1, t(Sk) = 0. Go to step 2. 



1. MaxF: descending number of followers, \Fj*\; 

2. MaxIF: descending number of immediate followers, \Fi\; 

3. MaxTime: descending task times, f,; 

4. MaxPW: descending positional weights, pwi — ti + Y^heF* OiJ 

5. MaxTimeL: descending task time divided by latest station, f,/L,-; 

6. MaxTimeSlack: descending task time divided by slack, f,/i, 

where £, and L, are the earliest and latest workstations (defined in Section |3.2.3[ ) to 
which a task can be assigned and still obey the preceden ce constraints and the desired 
cycle time c, and where Sj = L, — £, is the slack of task / 



With these priority rules, an algorithm for the SALBP-2 can be developed by 
searching for the smallest feasible c. SchoU and VoB ( 1996| l mention several search 



procedures such as the lower bound method and binary search. In this paper, we use 
the lower bound method in all our heuristics, since it finds the smallest possible c for a 
given SALBP-1 heuristic. 



3.2 Extension to the ALWABP 

The extension of the SALBP priority rules to the ALWABP faces two main difficulties. 
First, most of the priority rules rely on task execution times, which are not well defined 
for the case of the ALWABP since they depend on the designated worker. Second, there 
is no strategy to select the worker to be assigned to each workstation. In the following, 
we propose modifications to the original algorithm in order to cope with these two 
difficulties. 



3.2.1 Task priority rules 



From the six priority rules suggested by SchoU and VoB ( 1996| l, the first two can be used 
unchanged in ALWABP, because they do not depend on task execution times. Rules 
MaxTime and MaxPW do depend on task execution times but can still be adapted if one 
considers, for instance, minimum, maximum or averaged values of these times over all 
workers. In particular, for the original rule MaxTime based on task execution times, we 
considered ascending and descending orderings for each possibility. Finally, the last 



If £, = Li, Sj is set to some constant considerably smaller than one i SchoU and VoB 



1996 
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two rules depend on the earliest and latest stations to which a task can be assigned; we 
do decided to not apply these rules since the corresponding station bounds are rather 
weak for ALWABP. 

The proposed adaptions lead to eleven rules for the ALWABP derived from existing 
rules for the SALBP. We further propose five new rules that try to reflect the structure 
of the ALWABP by using the original task execution times per worker. The idea is 
to prioritize tasks when workers that quickly execute these tasks are being considered. 
The full set of priority rules are listed and described in the following, where w, S 
zrgmin^^,^^, t-tvi is some fastest worker in the execution of task / and = min„^wtwi, 
= max„^]vt„i, ti = LvvGW^vvi/l^l ^6 the minimum, maximum and mean execution 
time of task /, respectively. If task / is infeasible for worker w, i.e., t^t = °°, we set t^i = c 
for the purposes of computing and f;, where c is cycle time under consideration (see 
Algorithm |2]l. 

1. MaxF: descending number of followers, \F*\; 

2. MaxIF: descending number of immediate followers, |/^|; 

3. MaxTime^: descending minimum task times, f^^; 

4. MaxTime^: descending maximum task times, f,^; 

5. MaxTime: descending average task times. If, 

6. MinTime^: ascending minimum task time, f,^ ; 

7. MinTime+: ascending maximum task time, f,^; 

8. MinTime: ascending average task time, f,; 

9. MaxPW^: descending minimum positional weights, pwy = + Y.heF* '/T' 
10. MaxPW+: descending maximum positional weights, pwf = + Y.heF* ^t'^ 
IL MaxPW: descending average positional weights, pwi ~ ti+Y.heF* th', 

12. MinD(w): ascending difference to best worker, t„i — ty^.f, 

13. MinR(w): ascending ratio to best worker, t„i/t^,.i; 

14. MaxFTime(w): descending number of followers per time, \Fi\/t„i\ 

15. MaxIFTime(w): descending number of immediate followers per time, \F*\/t„i\ 

16. MinRank(w): ascending rank of worker's execution time, | {W G W \ t„ii < t„i] \ . 

Rules MinD(w) and MinR(w) prioritize tasks for which the current worker w is 
faster than other workers while MaxFTime(w) and MaxIFTime(w) balance two desired 
characteristics of a task: it should be executed quickly by the current worker w and it 
should make the largest number of tasks available. For the current worker w, rule 
MinRank(w) gives preference to a task with a short execution time compared to other 
workers, but does not depend on concrete execution times. The other rules follow the 
original rationale. 
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3.2.2 Worker priority rules 



Many of the task priority rules proposed above analyze the efficiency of a given worker 
for the task in comparison to the other workers. This suggests that Algorithm [T] needs 
to be modified to include a worker selection loop in which all possible workers are 
tested and, according to a worker selection criterion, one is chosen for the current 
workstation. This modification is reflected in Algorithm|2] We also add an outer while 
loop on cycle time tentative values, in order to present a self-contained algorithm for 
flie ALWABP-2. 



Algorithm 2 Station-oriented assignment procedure for the ALWABP-2 

1 : c = lower bound on the cycle time. 

2: while solution not found do 

3: Start at an empty workstation k = 0, f(5^.) = 0. 

4: Set the unassigned workers to Uw = W, set the unassigned tasks to t/ = A'. 

5: for each worker w in t/it. do 

6: Compute T,,,, the set of tasks that would be assigned to worker w according to the task 

priority rule used. 

7: end for 

8: Assign the worker w with the highest priority and its tasks T,,. to the station k. = 
U„\{w},U = U\T„,k = k+\. 

9: if U„ ^ then 

10: Go to step 4. 

11: else 

12: if f/ = then 

13: return c 

14: else 

15: c = c+l 

16: end if 

17: end if 

18: end while 



Algorithm |2] tests, for each workstation, all unassigned workers. In step 6, the 
algorithm computes the tasks that can be assigned to the current station if the worker 
under consideration is also assigned to that station. Note that these tasks depend on the 
priority rule, on the worker and on the tentative cycle time being considered. In step 
8, according to the selected tasks, it chooses the worker which will be assigned to the 
current workstation. Let ^ {W) = nm\wew twi be the minimum execution time of task ; 
among workers W. We propose three worker selection rules: 

• MaxTasks: descending number of assignable tasks, Ir^l; 

• MinBWA: ascending best worker assignment cycle MinBWA(f/ \ T„,Un, \ {w}) 
(Algorithm[3]). 

• MinRLB: ascending restricted lower bound, 

Lieu\T„ tr{U„\{w}) 
\U„\-l 
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Algorithm 3 Best worker assignment cycle (MinBWA) 



1: input: a set of tasks T and workers W 
2: Set = for all w in VK 
3: for each task ; in T do 

4: Determine the set of workers executing / fastest, W' := argmin,^,g)y 
5: Select a worker w' e W' of minimum load f^,/ = min,|,gvt" 

7: end for 

8: return the maximum cycle time c = max^g^y 



The first rule gives preference to the worker that was able to execute the largest 
number of tasks. The last two criteria compute estimates on the cycle time for the 
problem considering the yet unassigned workers and tasks, and give preference to the 
worker that reduces this estimate. They differ on the way the estimates are obtained. 
Rule MinBWA computes the cycle time that would be obtained if each remaining 
worker were assigned the tasks which he executes fastest (in case of ties, the task is at- 
tributed to the cuiTently less charged worker). Rule MinRLB computes a lower bound 
for the cycle time on the remaining workstations by dividing the sum of the minimum 
execution times of the remaining tasks by the number of remaining workers. 

Each of these three worker priority rules can be combined with one of the above 16 
task priority rules. Moreover, if one considers that the station-oriented assignment can 
be run in a forward (considering the original precedence graph) or backward (consid- 
ering the reversed precedence graph) manner, a total of 96 possible combinations are 
possible. The efficiency of these combinations is compared in Section|5] 



3.2.3 Bounds on the minimum cycle time 

In order to reduce the number of calls made to Algorithm|2]and, therefore, improve the 
efficiency of the method, one might investigate good lower bounds for the ALWABP- 
2. Valid lower bounds for the SALBP-2 can be used by substituting the task execution 
time in the original bound by the minimum task execution time over all workers f ,^ . 
With this strategy, we are able to use lower bounds LCI, LC2, and LC3 from[SchoU] 
and Becker, (2006) in our search procedures. LCI and LC2 are defined as 



LCI := max< max?,- 



ieN J 



LC2 := max( ^ f,7„^i_. | 1 < ^ < [(« - l)/mj 



The bound LC3 is obtained by destructive improvement ( Klein and Scholl \999) 
using bounds on the earliest station Ei{c) and latest station L,(c) for each task /, and 
is defined as the smallest c greater than or equal to an initial lower bound c such that 
Ei{c) < Li{c) for all tasks /. When using minimum task times, the station bounds are 
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defined as 



















L,(c) := m + 1 — 















For a given cycle time c, these bounds can be improved by preprocessing the 
ALWABP-2 instance as follows. For tasks i,k let 0{i,k) = {F*r]P^) U iF^*nP*) be 
the tasks that must be executed after / and before k or vice-versa, and let R{i,k) = 
0{i,k) U {i,k}. If some task / can only be executed by a given worker w, then another 
task k can be executed by the same worker w only if the total execution time of all tasks 
between / and k does not exceed c, i.e., Y.jeR(i.k) Uvj < c. Otherwise, task k is infeasible 
for worker w. In this latter case, a reduced instance for a cycle time c can be obtained 
by setting the execution time t^-k of all such tasks to oa^ making them infeasible for 
worker w. 

As a fourth bound we use the linear relaxation of model ([r|-([TO]i. For SALBP-2 
the linear relaxation of the ILP-model is known to be weak, but this bound proves to be 
useful for ALWABP-2, because the use of minimum execution times weakens the lower 
bounds LCI, LC2, and LC3. Since this bound is the most expensive computationally, 
it is used only in the genetic algorithm presented in the next section. 



4 A hybrid genetic algorithm for ALWABP-2 

In this section we propose a hybrid biased random-key genetic algorithm (HGA) based 
on the constructive heuristic of Section [3] Random-key genetic algorithms are par- 
ticularly useful for sequencing problems, and have been proposed by |Gongalves and| 
|de Almeida] ( |2004| l for the SALBP-1 and, by [Mendes et al| ( |2009l l for resource con- 
strained project scheduling. 



A genetic algorithm ( Goldberg 1989 Holland 1975 i is a population-based meta 



heuristic. Each individual in the population has a chromosome which codifies a solution 
to a problem instance. A chromosome consists of a collection of genes each of which 
can take a value among several alleles. A population generates offspring by crossover 
between two individuals and mutation. The probability of an individual to participate in 
a crossover is proportional to its fitness, making it more likely that fitter individuals pass 
genetic material to the next generation. A selection rule determines which individuals 
from the current generation and the offspring form the next generation. 

|Bean| ( |1994 ) introduced random-key genetic algorithms (RKGA). In a RKGA the 
chromosome is a vector in M" (where n is a problem-dependent parameter) and a de- 
coder maps it to a feasible solution. For example, in problems of sequencing jobs, 
the order of the genes' values defines a permutation of the jobs. The advantage of an 
RKGA is that the decoding process guarantees that each chromosome corresponds to a 
feasible solution and that the search can be done in a problem-independent way in the 
space W. 
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The evolutionary dynamics of a RKGA is as follows. All p individuals in the 
current population are sorted by their fitness value, and the best pe elite individuals 
are copied into the next generation. A RKGA replaces mutation by immigration of a 
small percentage p^ of random individuals into the next generation. The remaining p — 
Pe — Pr individuals are offspring. To generate offspring, two individuals are selected 
at random from the current population and are combined by uniform crossover. In 
uniform crossover, the offspring receives each allele independently with probability 
^ > 0.5 from the first parent and probability l—q from the second parent. 

A biased random-key genetic algorithm differs from a RKGA in the way it selects 
parents for crossover. The first parent is drawn randomly from the elite set, and the 
second one randomly from the remaining individuals. In this way, the offspring inherits 
with a higher probability an allele from an elite parent. A hybrid genetic algorithm 
applies a local search method to its individuals to improve their fitness. 

To construct a HGA for the ALWABP-2, we substitute the task priority rules pre- 
sented in Section 3 by explicit priorities. Instead of using a rule to define a priority 
pr(w, i) for each worker w and task /, the chromosome of an individual in the HGA is a 
matrix pm of priorities in [0, 1]. The fitness of an individual is the result of Algorithm 
2 run with these priorities. Algorithm|2]uses the bounds described in Section 3.2.3 



eluding the linear relaxation lower bound, when searching for a feasible cycle time. To 
find a feasible solution, it tries to allocate all tasks in a forward as well as a backward 
manner 

After decoding, a local search is applied to improve the solution. The local search 
uses two types of well-known moves ( |Scholl and Becker 2006| ): a shift of a task from 
one station to another, and a swap of two tasks between different stations. We allow 
also a sequence of two shifts, where the first one does not improve the solution. A third 
move type applied in the local search is a swap of two workers between two stations. 
All moves are applied only when they produce a feasible solution and are able to reduce 
the number of stations with load equal to the cycle time. 

The fitness of an individual is a pair (c, I), where c is the smallest cycle time found 
after local search, and I is the normalized total station time, i.e., the sum of all task 
execution times divided by mc. Fitness is ordered lexicographically. The second com- 
ponent allows the HGA to reduce, for a constant cycle time c, the station times of a 
solution. 

The initial population is seeded with individuals whose priorities are set according 
to the 16 proposed task priority rules. If the population size is less than 16 we select 
the p best individuals. For populations larger than 16, the remaining individuals of 
the population are initiaUzed randomly. The pseudo-code for the HGA is shown in 
Algorithm]?] 



5 Computational study 

In this section, we present a computational study of the proposed lower bounds, the 
constructive heuristics, and the hybrid genetic algorithm. In all experiments, we have 
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Algorithm 4 HGA for the ALWABP-2 



1: input: population sizes p, pe, Pr, crossover probability q; 

2: {We denote by /(/) the fitness of individual /} 

3: P:={(u,/(u)),...,(/i6,/(/i6))} where^e [0,1]I"'I^I'^I, 1 < A: < 16 is set using 
task priority rule k; 

4: P := PU{{in,f{iij)), {ip,f(ip))} where each i/^, 11 < k < p is a random ma- 
trix in [0, 

5: sort P by increasing cycle times f{ik) and keep the p best individuals; 

6: while stopping criterion is not satisfied do 

7: P' := . . . , {ipe,f{ipg)}', { elite passes to the next generation } 

8: tor k — I,. . . ,p~ pe~ Pr do 

9: select a random individual / from {i'l , . . . , ip^}; 
10: select a random individual j from {ip^+i ,.. .,ip}; 
11: o := crossover(/, j,^); 
12: localsearch(o); 
13: P':=P'U{io,f{o))}; 
14: end for 

15: for = 1, . . .,pr do 

16: P' := P'U {(/,/(/))} where / is a random matrix in [0, Ijl^l^l^l; 
17: end for 

18: P=P'; 

19: sort P by increasing cycle times 

20: end while 

21: return best individual u; 



used the instances available in the literature ( Chaves et al 2007 p They are grouped in 



four families: Roszieg, Heskia, Wee-Mag and Tonge, each one containing 80 instances. 
The instances were generated from the corresponding SALBP instances such that they 
contain 10 instances for each combination of five experimental factors at a low and 
a high level. The factors are the number of tasks, the number of workers, the order 
strengtlj^ the variability of the task execution time, and the number of infeasible task- 
worker pairs. For instances with low variability, the task execution times are drawn 
uniformly from the interval [l,f,], where f, is the task execution time as defined by the 
SALBP instance. When the variability is high, this interval was [l,3f,]. The number 
of infeasible task-worker pairs is 10% and 20% on the low and the high level, respec- 
tively. The main characteristics for each group of instances (number of tasks, number 
of workers and the order strength of the precedence network) are listed in Table [T] 

Our results are given as relative deviations over the objective value of the optimal 
solution or the best known value. For instance families Heskia and Roszieg the optimal 



solutions are known (for a complete table see Blum and Miralles (201 1 1). For instance 



families Wee-Mag and Tonge we compare with the best known value defined as the 



^The instances can be found at http : //www. f eg .unesp .br/'chaves/Arquivos/ALWABP . zip. 
'The order strength is defined as the fraction of present precedence relations compared to the maximum 
possible, i.e., 2\E*\/ (ii{n — 1)). 
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Table 1: Instance characteristics. 



Family 




\w\ 


Order Strength 


BKV/LB (%) 


Roszieg 


25 


4 (groups 1-4) or 6 (groups 5-8) 


71.67 


16.75 


Heskia 


28 


4 (groups 1-4) or 7 (groups 5-8) 


22.49 


57.91 


Tonge 


70 


10 (groups 1-4) or 17 (groups 5-8) 


59.42 


57.91 


Wee-Mag 


75 


11 (groups 1-4) or 19 (groups 5-8) 


22.67 


76.70 



minimum of the best values found by a clustering search, a tabu search and an iterated 



beam search as published in the literature (Chaves et al 2009 Moreira and Costa 2009 



[Blum and Miralles 201 1) . All computation times reported are in seconds of real time. 



To be able to estimate the quality of the lower bounds and the best known values, we 
report in the last column of Table[T]the average relative deviation of the optimal solution 
or best known value from the lower bound. In comparison over all 320 instances lower 
bound LCI was maximum in 118 cases, LC2 in none, LC3 in 162 cases and the linear 
relaxation of the ILP model ([T])-([TOl) in 235 cases. 

5.1 Constructive heuristic methods for ALWABP-2 

For the SALBP, it is known that breaking ties between tasks of the same priority can 



have a significant effect on the solution quality (Talbot et al 1986 1. Therefore we tested 
several tie-breaking rules for both task and worker selection. For selecting tasks, we 
used as a first-level tie-breaker the descending number of immediate followers (MaxIF) 
and as a second-level tie-breaker the task execution time for the current worker (f».,). 
For the worker selection, as a first-level tie-breaker, we used MinRLB for rules Max- 
Task and MinBWA and MaxTask for rule MinRLB. As a second-level tie-breaker, we 
used the workstation idle time. Both for tasks and workers, the original index of the 
task or worker was considered as a third-level tie-breaker, in order to make the heuris- 
tics deterministic. 

Another important implementation aspect is the fact that some of the priority rules 
can be applied statically or dynamically with respect to the worker selection loop. For 
instance, task priority rules MinD(w) and MinR(w) as well as all task priority rules that 
rely on the calculation of a minimum, maximum or average value of a parameter over 
all workers are prone to changes as the set of available workers is reduced. Preliminary 
tests showed that the dynamic strategy using only the available workers presented better 
results and was used in all subsequent tests. For the sake of simplicity and due to the 
computational efficiency of Algorithm 2, we used LCI to obtain the lower bound on 
the cycle time in line 1 of Algorithm 2. 

Under the considerations above, all combinations of task and worker priority rules 
were considered. The heuristics were implemented in C, using the compiler gcc 4.4, 
under the Linux Ubuntu operating system. For the tests we used a PC with a Core 2 
Duo 2.2 GHz processor and 3 GB of main memory. 

The first conclusion that could be drawn was that for all possible task selection 
rules, worker selection rule MinRLB consistently obtained better results when com- 
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pared to rules MaxTasks and MinBWA, which were, therefore, discarded. The results 
obtained for the 16 task priority rules, worker selection rule MinRLB and strategies 
backward and forward are presented in Table |2] In the table, we present the average 
and maximum deviation of the heuristic solution with respect to the best known values 
and the average and maximum computation time. 

From the results in Table |2] we can conclude that the heuristics can obtain feasi- 
ble solutions for the ALWABP in very small computation times (not more than 0.27 
seconds in the worst case). The quality of the obtained solution is sensitive to the task 
priority rule used, with the rules based on the positional weights presenting the best re- 
sults. Indeed, criterion MaxPW^ yielded average deviations of 17.0 percent and 26.2 
percent for the forward and backward searches, respectively. The diversity of the task 
priority rules can be evaluated by considering, for each instance, the best obtained so- 
lution over all rules. This computation yields an average deviation of 9.6 percent in 
less than one second of average computation time. The application of the local search 
described in Section[4]presented only marginal gains. 

These results are competitive with those obtained by the clustering search algorithm 
of Chaves et al (2009i which obtained average deviation values of 21.1 percent and 
with the tabu search of Moreira and Costa] ( |2009| l which obtained average results of 
26.8 percent, in larger computation times. Indeed, both meta-heuristic methods can 
benefit from an initial solution such as those obtained with our proposed strategy. We 
ran both methods initialized with the solutions obtained with criterion MaxPW^ and 
forward allocation. The clustering search algorithm and the tabu search algorithm then 
presented average deviations of 1 1 .7 percent and 9.9 percent respectively. These results 
can be further improved if the best solution (over all priority criteria) is used as initial 
solution. In this case, the clustering search and the tabu search yielded average results 
of 7.8 percent and 7.4 percent, respectively (see Table [3]l. These results suggest the 
advantage of quickly obtaining good initial solutions such as those presented in this 
article. 



5.2 A HGA for ALWABP-2 

Based on the results of the previous section we chose MinRLB as the worker priority 
rule for the HGA. In a preliminary experiment on three selected instances of all four 
families, we tested the performance of the HGA on populations of size P = 40, 70, 100, 
and for uniform crossover probabilities of p ~ 0.5, 0.6, 0.7, 0.8, 0.9. Based on this 
experiment, we chose a population size of P = 100 and a crossover probability of 
p = 0.5. We ran the HGA with these parameters 20 times with different random seeds 
on all 320 instances. The execution stopped after 200 iterations, or 100 iterations 
without improvement of the incumbent. We used a PC with a 2.8 GHz Intel Core 
17 930 processor with 3 GB of main memory for the experiments. 



Tableplcompares the results of the HGA with the clustering search (CS) of Chaves 
let al| ( |2009| , the tabu search (TS) of [Moreira and Costa] ( |2009| l, and the iterated beam 
search (IBS) of ]Blum and Mrralles '(201 1). For each group of the four instance families, 
and for each of the four meta-heuristics, we report the relative deviation of the best 
value found from the best known value, the relative deviation of the average value from 
the best known value, the average total computation time, and the average computation 
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Figure 1 : Main effects (upper row) and most significant double interactions (lower row) 
for HGA and IBS of the five experimental factors on average relative deviations. 



time to find the best solution. We report only the solution value for TS, since the 



method is deterministic, and only the time to find the best value for IBS, since |Blum 



and Miralles (2011 1 do not report total computation time. The results of CS and of 



the TS have been obtained using the best constructive heuristic to generate the initial 
solution. 

All meta-heuri sties have similar computation times. The results of the CS have 
been obtained on a PC with a Pentium 4 2.6 GHz with 1 GB main memory, the TS 
results on a PC with a 1.66 GHz Intel Core 2 Duo T5450 with 3 GB main memory, 
and the results of the IBS on a PC with a 2.6 GHz Pentium 4 processor and 1 GB of 
main memory. The computation time on these architectures may vary about a factor of 
two, but we consider the total computation time of all methods to be reasonable for an 
NP-hard optimization problem. 

Concerning the solution quality, for the smaller instance families Heskia and Roszieg, 
all methods obtain good results, with a relative deviation of the average value of at most 
1.5 percent (CS) and a relative deviation of the best value of at most 0.4 percent (CS). 
IBS and TS solve these instances optimally. The small relative deviations of the HGA 
come from two of the 80 instances, which were not solved to optimality. On the larger 
instance families Wee-Mag and Tonge, IBS and the HGA always outperform CS and 
TS, with average relative deviations in each group of instances at least four percent 
better The HGA obtains better best values in eight of the 16 groups corresponding to 
73 of the 160 instances. The HGA is able to improve the average values in eight of the 
16 groups. Table |4] summarizes the number of instances in which IBS and HGA obtain 
a better average and best value. 

A statistical analysis shows that the most significant factors influencing the average 
solution quality of HGA and IBS are the number of tasks and workers, and the order 
strength. Figure [T] shows the main effects and interaction plots for the two most sig- 
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Table 4: Comparison of IBS and HGA. 





IBS better 


Ties 


HGA better 


Avg. 


93 


143 


84 


Best 


48 


199 


73 



nificant double interactions for HGA and IBS. The solution quality of both methods 

strongly depends on the number of tasks. IBS performs better for low order strength, 
and the HGA performs significantly better when the number of workers is low. These 
dependencies are more pronounced for a high number of tasks. 

In summary, the approach using a HGA seems adequate for obtaining good solu- 
tions, with results comparable to or better than the existing heuristics. In particular it 
obtains significantiy better solutions for a low number of workers. Being a relatively 
simple method, optimizing the task priorities of a station-based assignment scheme in- 
side a generic hybrid genetic algorithm, the HGA may be an interesting alternative for 
solving the ALWABP-2. 

6 Conclusions 

In this paper, we have proposed a constructive heuristic framework based on task and 
worker assignment priority rules for the assembly line worker assignment and balanc- 
ing problem. In a series of computational tests, the approach proved to be fast and 
efficient both when evaluated as a stand-alone method as well as when the obtained 
solutions were used to improve the convergence of more elaborate meta-heuristics. 
Moreover, the strategy was used as a solution decoder within a hybrid genetic algo- 
rithm which was also proposed and tested, obtaining results that were comparable to 
the best known methods available in the Uterature. 
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